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Abstract 

A  model  which  incorporates  the  influence  of 
electrode  surface  conditions,  gas  pressure,  and 
charging  rate  on  the  voltage  stability  of  high  energy 
spark  gaps  is  discussed.  Implications  of  the  model 
include  changes  in  the  width  of  the  self-breakdown 
voltage  probability  density  function  as  the  secondary 
emission  characteristics  of  the  cathode  are  modified 
by,  for  example,  oxide  and  nitride  coatings  and/or 
deposits  from  the  insulator.  The  model  indicates  that 
a  narrow  self-breakdown  voltage  distribution  requires 
a  source  of  electrons  near  the  cathode  surface,  which 
could  be  provided  by  UV  photo-illumination  of  the 
cathode.  In  addition,  the  model  provides  an  extremely 
useful,  and  physically  reasonable  framework,  from 
which  the  properties  of  spark  gaps  under  a  wide 
variety  of  experimental  conditions  may  be  evaluated. 

Both  experimental  and  theoretical  results  are 
presented . 

Introduction  V 


Low- jitter,  triggered  spark  gaps  are  needed  for  a 
wide  variety  of  switching  applications  including 
fusion  machines  [1],  weapons  systems,  and  high  energy 
physics  experiments.  To  achieve  low  jitter,  the 
switch  should  be  triggered  as  close  to  the  self-break- 
down  voltage  as  possible;  thus,  an  Ideal  switch  should 
have  a  delta  function  for  the  self-breakdown  voltage 
probability  density  function,  In  actual  operation  the 
self-breakdown  voltage  will  be  somewhat  erratic  and  in 
most  cases  "prefires",  or  breakdown  voltages  which  are 
significantly  less  than  the  mean,  will  occur.  The 
self -breakdown  voltage  density  functions  and  the 
respective  distributions  for  these  cases  are  shown  in 
Figure  1 . 

Numerous  studies  [2-5]  have  shown  that  the  choice 
of  gas,  electrode  and  insulator  material  can 
significantly  influence  the  width  and  shape  of  the 
actual  voltage  density  function.  More  specifically, 
several  studies  [6-8]  have  suggested  a  correlation  of 
the  statistical  distribution  in  the  self  breakdown 
voltage  of  a  spark  gap  and  the  properties  of  the 
cathode  surface,  including  its  microstructure.  The 
data  have  been  interpreted  in  terms  of  models  that 
consider ! 


Fig.  1.  The  Self-Breakdown  Voltage  Probability 

Density  Function,  py(v),  and  its  Distribution 
Function,  Fy(v)  for  a)  Ideal  Spark  Gap,  b) 
Actual  Spark  Gap,  c)  Actual  Spark  Gap  with 
Prefires . 

has  been  satisfied.  The  theoretical  model  and  the 
experimental  results  presented  here  Includes  both  of 
these  mechanisms  by  which  the  cathode  surface  can 
affect  the  statistical  distribution  in  the  breakdown 
voltage  and,  in  addition,  includes  the  field 
enhancement  effects  on  the  cathode  surface  in  a  new 
way. 


1  )  the  effect  of  the  cathode  surface  properties 
on  the  primary  election  emission  current,  ie,  at  the 
cathode  [6,9,10],  and 

2)  the  effects  of  the  field  enhancement  due  to 
cathode  surface  microstructure  on  Townsend's  first 
ionization  coefficient,  [6,11,12]. 

Both  models  include  the  concept  of  waiting-for-an- 
electron,  in  that  breakdown  is  assumed  to  occur  when 
the  first  electron  appears  at  the  cathode  after  a 
breakdown  condition  (usually  the  Townsend  condition) 
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Theoretical  Model 

General  Case 

Consider  a  spark  gap  subjected  to  a  monotonically 
increasing  applied  voltage,  v(t).  Denote  the 
breakdown  voltage,  a  random  variable,  as  V.  The  field 
enhancement  factor,  M,  on  the  cathode  surface  is  also 
considered  to  be  a  random  variable.  (The  underlying 
sample  space  is  the  geometrical  surface  of  the 
cathode. )  The  random  variable  M  is  characterized  by  a 
probability  density  function,  p^(m).  A  basic 
assumption  of  the  model  is  that  the  gap  breaks  down 
when  an  electron  is  born  at  a  site  on  the  cathode 
surface  where  M  is  as  large  as  or  larger  than  the 
value  that  satisfies  the  breakdown  condition  (Townsend 
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perhaps)  at  the  particular  voltage  applied.  We  denote, 
this  threshold  value  of  the  field  enhancement  as 
jnt  ( v ) .  Physically  we  expect  that  m^fv)  is  a 
monotonically  decreasing  function  of  v  Omt(v)/3v<0)  , 
an  increasing  function  of  pressure  (3nq-(v)/3p>0) ,  and 
that  mj-(0)=  co  and  m^-(o3)=0.  The  actual  calculation  of 
mt(v)  has  been  derived  and  explained  elsewhere  [8,13], 

We  now  calculate  the  probability,  p^,  that  the 
gap  breaks  down  during  the  time  between  t  and  t  +  At: 
and  hence  at  a  voltage  between  v  and  v  +  Av.  For  At 
small,  the  probability  that  an  electron  is  born 
between  t  and  t  +  At  at  a  site  where  M  takes  a  value 
between  m  and  m  +  Am  is 


ie(m,v(t) ) 

-  At  p^nOAm  (1) 

e 

The  quantity  e  is  the  charge  on  an  electron  and  ie, 
the  primary  electron  current  generated  at  the  cathode, 
is  allowed  to  depend  on  both  the  field  enhancement,  M, 
at  the  cathode  surface  and  on  the  applied  voltage, 
v(t).  Thus,  the  probability  that  the  gap  breaks  down, 
between  t  and  t  +  At  is 
At  03 

pt(At)  «  —  J  ie(m,v(t)  )p[q(m)dm  (2) 

6  mt(v) 


It  can  be  shown  [13]  that  this  leads  to  the  expression 
for  the  voltage  density  function  py(v) : 


Pv(v) 


M  v) 


-  exp 

v' 


\(ri)dri 

v' 


where 


X(v)  =  -  J  ie(m,v)p^(m)dm 
6  mt(v) 


(3) 


(4) 


Special  Cases 


To  proceed  further,  consider  two  special  cases  of 
the  model.  First,  suppose  that  ie(m,v)  is  constant 
so  that  field  enhancement  distribution  effects  from 
the  cathode  surface  microstructure  and  waiting-f or-an- 
electron  effects  are  the  primary  physical  mechanisms 
included  in  the  model.  This  circumstance  is  likely  to 
hold,  for  example,  when  the  cathode  is  illuminated 
with  sufficiently  intense  ultraviolet  radiation  so 
that  any  field  emission  current  is  dwarfed  by 
photoelectric  current,  which  should  be  independent  of 
M  and  v.  If  ie  =  ieo ,  a  constant,  then  Eq.  (6)  gives 


Fjq(mt(v) ) 


e  v'py(v) 
ieo  l~Fy(v) 


(7) 


where  F{q(m)  is  the  cummulative  probability  distributi¬ 
on  for  the  random  variable,  M.  If  we  know  m^-(v),  we 
can  determine  Fjq(m)  by  plotting  Fjqtm^fv))  vs  irq-(v). 
For  this  special  case,  therefore,  it  is  possible,  in 
principle,  to  deduce  Fjq(m)  from  py(v)  (self  breakdown 
voltage  histogram)  under  a  given  set  of  conditions  and 
thus  predict  p\/(v)  (or  Fy(v))  for  a  different  v'  or 
gas  pressure  (which  affects  irq-fv)),  for  example.  For 
this  special  case,  Eq.  (5)  becomes 


Fy(v)  *  1-exp 


e 


v 

J 


o 


[  1 — ( mt ( rj )  )  ]dn 

v’ 


(8) 


Consider  now  a  second  special  case  in  which  pjq(m) 
=  S(tn  -  m0 ) ,  where  {(•)  is  the  Dirac  delta  function 
and  m0  is  a  constant.  In  this  case  the  field 
enhancement  is  assumed  to  be  uniform  (that  is, 
sufficiently  characterized  by  its  mean  value  rather 
than  its  distribution)  so  that  the  primary  effects 
included  are  the  voltage  dependence  of  the  primary 
electron  current,  ie,  and  waiting-f or-an-electron.  In 
this  case  Eq .  (5)  becomes 


It  is  easy  to  see  that 


v 

Fv(v)  =  /  pyUMS 
o 


X  ( n )  dn 

v  ( n ) 


(5) 


where  Fy(v)  is  the  cumulative  probability  distribu¬ 
tion  for  the  random  variable,  V.  Note  that  v'is  to  be 
evaluated  at  v,  and  hence  can  be  considered  as  a 
function  of  v(v'=v'(v)). 


Fy(v)  =  1-exp 


1  v  i@(m0 ,q)dn 


(9) 


where  v^-  is  the  threshold  voltage,  and  m^  ( v-t- )  =m0  , 
while  Eq.  (6)  becomes 


ie(m0 , v) 


ev'py(v) 

1-Fv(v) 


(10) 


If  v(t)  is  a  ramp,  then  v'=v0',  a  constant.  If 
v(t)  is  an  RC  charging  waveform,  then  v'=(v0-v)/RC 
where  v0  is  the  charging  voltage. 

These  equations  for  py(v)  and  Fy(v)  represent  the 
model's  prediction  of  the  voltage  self-breakdown 
statistics  for  the  gap  in  terms  of  the  primary 
electron  current,  ie(m,v);  the  probability  density 
function  p^(m)  of  the  cathode  field  enhancement, 
considered  as  a  random  variable;  the  threshold  field 
enhancement,  mt(v),  expressed  as  a  function  of  the 
voltage  applied  in  the  gap;  and  the  time  derivative, 
v',  of  the  applied  voltage.  From  these  equations  it 
is  easy  to  show  that 


1  m  v'py(v) 

-  J  ie(m, v)pM(m)dm  “  -  (6) 

e  _  ,  .  1  -  Ft,! \j) 


Notice  that  the  right  hand  side  of  Eq .  (6)  can  be 
measured  experimentally. 


Experimental  Arrangement 

The  experimental  arrangement  and  the  system 
diagnostics  used  to  test  the  theoretical  results  are 
shown  in  Figures  2  and  3.  The  construction  of  this 
facility  and  the  development  of  the  modeling  software 
is  described  elsewhere  [8,14],  The  test  circuit  shown 
in  Figure  3  consists  of  a  high  energy  ( 2kJ )  pulse 
forming  network  (PFN)  and  a  low  energy  (<lmJ)  RC 
probing  circuit.  The  PFN  is  used  to  generate  an 
electrode  surface  which  is  characteristic  of  a  high 
energy  switch.  The  RC  probing  circuit  is  used  to 
generate  the  voltage  distributions  with  a  low  energy, 
low  current  pulse  so  that  the  steady  state  temperature 
is  reached  prior  to  each  shot.  This  low  energy 
circuit  is  also  used  so  that  the  surface  microstruc- 
ture  will  not  be  altered  significantly  from  shot  to 
shot.  The  pressure  in  the  spark  gap  could  be  raised 
to  3.5  atmospheres  and  the  voltage  ramp  rate  could  be 
varied  from  3  to  60  kV/sec  by  changing  Rc.  A  5  Watt  UV 
lamp  was  used  to  generate  additional  electrons  at  the 
cathode  surface  when  needed. 
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Fig.  2.  Experimental  Arrangement  and  System 
Diagnostics . 


Fig.  3.  Test  Circuits. 


A  testing  sequence  consisted  of  firing  2000-7000 
shots  at  high  energy,  waiting  for  approximately  1  hour 
for  the  electrode  to  cool,  and  proceeding  with  several 
series  of  500,  low  energy  shots  with  different  ie,  v' 
and  pressure.  A  typical  electrode  surface  generated 
by  the  high  energy  pulses  for  the  case  of  304 
stainless  steel  run  in  1  atmosphere  of  nitrogen  gas  at 
a  gap  separation  of  5mm  is  shown  in  Figure  4. 
Examination  of  the  electrode  surface  after  application 
of  the  low  energy  pulses  Indicates  that  no  significant 
changes  had  occurred  which  might  alter  the  breakdown 
statistics . 


Fig.  4.  Stainless  Steel  Electrode  Run  in  Nitrogen,  a) 
Cathode  Surface-Top  View,  b)  Cathode  Surface- 
Side  View. 


Experimental  Results 

Several  experiments  were  performed  to  verify  the 
nodel's  predictions  for  the  effect  of  ie,  v',  and 
presence  on  the  probability  density  function  py(v).  In 
the  first  experiment,  an  UV  source  was  used  to 
generate  a  continuous  supply  of  photoelectrons  at  the 
surface  of  a  stainless  steel  electrode  in  air.  Figure 
5  shows  that  without  UV,  the  density  function  is  very 
broad  indicating  that  the  cathode  surface  is  a  very 
poor  emitter  of  electrons.  However,  with  the  UV 
source  on,  the  density  function  collapsed  to  the 
lowest  value  of  breakdown  voltage.  This  result  is 


significant  for  at  least  two  reasons.  First,  it 
supports  the  waiting-f ox-an-electron  concept  as  one 
mechanism  responsible  for  statistical  variation  in  the 
self-breakdown  voltage,  and  second,  it  provides  an 
externally  controllable  experimental  "switch"  where 
the  effect  of  waiting-f or-an-electron  can  be  turned  on 
or  off.  The  behavior  observed  is  consistent  with  Eq . 
(5). 


Fig.  5.  Self  Breakdown  Voltage  Probability  Density 
Function  For  Stainless  Steel  In  Air  With  and 
Without  UV. 

A  second  experiment  consisted  of  varying  the 
voltage  ramp  rate  v'  from  3  kV/sec  ("slow"  ramp)  to  30 
kV/sec  ("fast"  ramp).  According  to  the  model  (Eq, 
(5)),  if  you  are  waiting  for  an  electron  to  appear, 
then  the  faster  the  ramp  rate,  the  higher  the 

breakdown  voltage  will  be  when  the  electron  appears 
and  thus  the  greater  the  scatter  in  the  density 

function  py(v).  Figure  6  shows  that  this  effect  was 
indeed  observed.  Also,  from  Eq .  (5)  the  density 

function  for  the  slower  ramp  rate  could  be  theoretic¬ 
ally  calculated  from  the  data  for  the  fast  ramp  rate. 
Figure  6  shows  this  result  for  the  assumption 

ie  =  ie0,  a  constant.  The  result  is  fair  at  best: 

which  indicates  that  for  better  agreement  a  more 
realistic  expression  for  ie,  perhaps  ie(m,v),  would 
have  to  be  used. 


Fig.  6.  Self  Breakdown  Voltage  Probability  Density 
Function  For  Different  Charging  Rates. 

Previous  work  [12,15-17]  has  shown  that  with  the 
presence  of  cathode  microstructure,  an  increase  in 
pressure  can  lead  to  significant  deviations  from  the 
Paschen  curve  breakdown  voltage  if  the  product  of  the 
protrusion  height  and  the  pressure  is  greater  than  a 
gas  dependent  threshold.  Avrutskil  (7)  stated  that  an 
increase  in  pressure  should  lead  to  an  increase  in 
scatter  in  the  breakdown  voltage  but  no  data  were 
given.  Thus,  in  order  to  understand  the  effect  of 
pressure  on  the  breakdown  voltage  statistics  for  a 
surface  with  large  protrusions,  a  brass  sample  was 
generated  and  the  breakdown  voltages  were  recorded  for 
pressures  up  to  3.5  atmospheres.  (Earlier  work  in 
electrode  erosion  had  shown  that  brass  electrodes  in 
high  energy  operation  could  form  protrusions  up  to 
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500  urn  (14).)  Figure  7  clearly  shows  an  increase  in 
scatter,  especially  at  the  low  end,  in  the  density 
function  py(v)  for  higher  pressures. 


Fig.  7.  Self  Breakdown  Voltage  Probability  Density 
Function  For  Different  Pressures. 

The  pressure  data  was  found  to  be  of  even  greater 
importance  for  analyzing  the  different  ie  cases  which 
were  studied  theoretically.  Figure  8  shows  theoret¬ 
ical  plots  of  the  quantity  v'py(v)/( l-Fy(v) ) ,  for  the 
three  physical  cases  discussed  earlier:  a)  ie  = 
ie(m0,v)  m0  is  a  constant  over  the  entire  surface  (Eq. 
(8)),  b)  ie  »  ieo  is  a  constant  (Eq.  (9))  and  c)  the 
most  general  case,  ie  =  ie(m,v)  which  assumes  a 
distribution  of  surface  field  enhancements  (Eq.  (5)). 


Fig.  8.  Theoretical  Plots  of  v'py(v)/( l-Fy(v) )  For  a) 
ie  “  ie(m0,v),  b)  ie  ■  ieo  and  c)  ie  » 
ie(m, v) . 

Case  a)  illustrates  that  if  no  distribution  exists  in 
M  (Eq.  (10))  then  an  increase  in  pressure  will 
correspond  simply  to  a  higher  emission  current  because 
of  the  higher  breakdown  voltage,  which  is  typical  for 
a  field  dependent  Schottky  or  Fowler-Nordheim  emission 
mechanism  [9J.  However,  in  case  b)  for  a  fixed 
voltage,  the  increase  in  pressure  has  the  effect  of 
raising  the  threshold  mt  required  for  breakdown  which 
raises  F^fm^lv))  and  thus  the  function  v'py(v)/ 
(l-Fy(v)),  is  multivalued  and  decreases  with 
increasing  pressure  (Eq.  (7)).  For  case  c)  v'py(v)/ 
(l-Fy(v)),  is  also  multivalued  and  decreases  with 


increasing  pressure  in  the  following  way(Eq.  (6)). 
For  a  fixed  voltage  and  assuming  the  surface  features 
do  not  change  with  pressure,  the  integrand  is  constant 
with  increasing  pressure.  However,  the  lower  limit  on 
the  integral,  namely  mt(v),  increases  with  increasing 
pressure  which  has  the  effect  of  reducing  the  value  of 
the  function. 

Figure  9  is  a  plot  of  the  function  v'py(v)/ 
(l-Fy(v)),  from  experimental  data  for  the  pressure 
data  of  Figure  7.  From  this  plot  it  is  clearly  seen 
that  the  experimental  data  are  inconsistent  with  the 
theoretical  results  for  case  a)  (a  constant  m 
surface).  Thus,  the  effect  of  a  distribution  in  field 
enhancements  must  be  considered  in  the  analysis  of  the 
breakdown  statistics. 


Fig.  9.  Experimental  Plots  of  v 1 py( v) /( 1-Fy( v) ) . 


Conclusion 


Theoretical,  and  experimental  results  show: 

1)  The  spread  in  self-breakdown  voltages  in  a 
spark  gap  is  a  function  of  the  charging  rate  (v' )  and 
the  cathode  surface  properties  which  determine  the 
electron  emission  current  ie  and  the  distribution  of 
field  enhancement  sites  F^(m). 

2)  Changing  ie  provides  a  practical  method  for- 
reducing  the  width  of  the  self -breakdown  voltage 
density  function.  This  can  be  accomplished  with  an 
external  UV  source  or  perhaps  with  an  electron 
emission  agent  introduced  into  the  cathode  material 
[18]. 

3)  At  sufficiently  high  pressures  the  cathode 
microstructure  has  adverse  effects  on  the  voltage 
density  function. 

4)  The  function  v'py(v)/( l-Fy(v) ) ,  which  can  be 
computed  directly  from  self-breakdown  voltage  data,  is 
useful  for  determining  the  nature  of  ie  for  a  given 
set  of  conditions. 

Some  questions  of  interest  which  remain  to  be 
answered  include: 

1)  Under  what  conditions  does  the  cathode  micro¬ 
structure  significantly  affect  the  voltage  stability 
of  a  spark  gap? 

2)  What  is  the  physical  mechanism  which  produces 
the  electron  emission  current  and  can  it  be  artific¬ 
ially  stimulated? 

3)  What  is  the  mechanism(s)  which  produces  pre¬ 
fires  and  can  it  be  controlled? 
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